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ABSTRACT 

We present the first mid-infrared Spiteer/Infrared Spectrograph (IRS) observations of powerful radio 
galaxies at z>2. These radio galaxies, 40-^23.56 (z=2.48) and 6CJ1908+7220 (z=3.53), both show 
strong mid-infrared continua, but with 6C J1908-I-7220 also showing strong PAH emission at rest- 
frame 6.2 and 7.7 /xm. In 4C -1-23.56 we see no obvious PAH features above the continuum. The 
PAH emission in 6C J1908-I-7220 is the amongst the most distant observed to date and implies that 
there is a large instantaneous star formation rate (SFR). This is consistent with the strong detection of 
6C J1908+7220 at far-IR and sub-mm wavelengths, indicative of large amounts of cold dust, '^10^ Mq. 
Powerful radio galaxies at lower redshifts tend to have weak or undetectable PAH features and typically 
have lower far-IR luminosities. In addition, 4C -1-23.56 shows moderate silicate absorption as seen in 
less luminous radio galaxies, indicating r9.7^m=0.3 ± 0.05. This feature is shifted out of the observed 
wavelength range for 6C J1908-I-7220. The correlation of strong PAH features with large amounts of 
cold dust, despite the presence of a powerful AGN, is in agreement with other recent results and implies 
that star formation at high redshift is, in some cases at least, associated with powerful, obscured AGN. 
Subject headings: galaxies: active — galaxies: high-redshift — galaxies: evolution - radio galaxies: 
individual(4C -h23.56, 6C J1908-h7220) 



1. INTRODUCTION 

There is now clear evidence that high redshift, pow- 
erful radio galaxies (HzRGs; z > 1 and Lsghz > 
10^^ WHz-i) have massive, ~ 5 x 10" M©, stellar hosts 
(Seymour et al. 2007; see also Lilly & Longair 1984; De 
Breuck et al. 2001; Rocca-Volmerange et al. 2004) and 
that the stellar mass of t hese hosts is related to the ra- 
dio p ower of the HzRGs (|Lacv et al.ll2000l : IWillott et all 
|2003[) . The strong radio and mid-IR continuum lumi- 
nosities of these rare sources imply hig h accretion rates 
onto their central massive black holes (|0gle et al.l 120061 : 
ISevmour et al.l l2007f) . There is also evidence for ex- 
tremely high SFRs from the strong sub-mm emission seen 
in some HzRGs ( lArchibal d et al.l l2001t iReuland et all 
l2004t iGreve et al.i r2006). as well as from the spectacu- 
lar (> 1 00 kpc) luminous Ly g haloes observed in many 
HzRGs (|Reuland et all 120031 : lyTllar-Martin et al.ll2003f ). 
As the time scales of both the radio jets and star for- 
mation are very short (< 0.1 Gyr), and as HzRGs are 
rare objects, it seems we are observing these sources at a 
special phase during a strong spurt of coeval black hole 
and galaxy growth, but when both processes are dust- 
enshrouded. 

Aside from sub-mm/far-IR observations, the only 

^ Spitzer Science Center, California Institute of Technology, 1200 
East California Boulevard, Pasadena, CA 91125, USA. [email: 
seymourOipac . caltech.edu] 

^ European Southern Observatory, Karl Schwarzschild Strai3e, 
D-85748 Garching, Germany. 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138, USA. 

* Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, CA 91109, USA. 

^ Astronomisches Institut, Ruhr-Universitat Bochum, Univer- 
sitatsstr. 150, Gebaude NA 7/173, D-44780 Bochum, Germany. 

^ Departments of Physics Sz Astronomy, Marian E. Koshland In- 
tegrated Natural Science Center, 370 Lancaster Avenue, Haverford 
College, Haverford, PA 19041, USA. 



other method to obtain an approximate measure of the 
SFR in HzRGs is from the mid-IR PAH emission fea- 
tures. This diagnostic is known to be inexact, especially 
in AGN, but UV/optical SFR diagnostics (e.g. narrow 
lines) are particularly difficult for HzRGs, as are radio/X- 
ray diagnostics. All these observables are known to be 
significantly contaminated by the powerful AGN. A cor- 
relation between PAH an d sub-mm luminosity has been 
found in local starbursts (jHaas et al.ll2002l : iPeeters et al.l 
120041) showing that PAH emission can be used as an ap- 
proximate tracer of star formation. There are now sev- 
eral results in the literature linking PAH emission and 
the presence of cold dust associated with star formation 
(Lacy et al. 2007; Pope et al. 2007; Sajina et al. 2008). 

This letter presents Spitzer/TRS spectra of two z > 

2 HzRGs from o ur sample of 69 HzRGs studied 

bv ISevmour et al.l (I2007D : 4C-h23.56 at z=2.48 and 
6C J1908-f-7220 (less commonly known as 4C -1-72.26) at 
z=3.53. Very few HzRGs have photometry at 24 /xm and 
even fewer are bright enough to be easily targeted by 
IRS. We select two of three HzRGs with S'24/im > 1 mJy 
which are also detected at 70 fim. While these HzRGs 
may not be representative of the entire population, even 
this small sample shows diverse mid-IR properties. 

We also fit the total IR, 1 - 1000 ^m, SEDs of these 
HzRGs and compare the model fits to the IRS spec- 
tra. We adopt the concordance ACDM cosmology with 
nM=0.3, and Hq^IO kms-^Mpc'^ 

2. THE SPITZER MID-IR SPECTRA 

The IRS spectra of 4C -1-23.56 were obtained as part of 
our Cycle 4 GTO program (PI G. Fazio) and consisted 
of 12 ramp cycles of 120 s each with just the LLl module 
(~ 20 - 36 /im). 6C J1908+7220 was observed as part 
of a Cycle 1 GO program (PI M. Haas) and the data 
include 3 cycles x60s with the SLl module, and 3 cycles 
xl20s with both the LLl and LL2 modules. We used 
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the interactive IRSCLEAN"'^ task to remove bad pixels in 
the Spitzer post-BCD images and extracted the spectra 
with SPICE^ using 'optimal' extraction. The spectra were 
smoothed by different amounts in each module varying 
from 1 fim at the shortest wavelengths to 0.3 /im at the 
longest wavelengths. 

Figuredlpresents the IRS spectra and Figure[2]presents 
full infrared spectral energy distributions (SEDs) of both 
HzRGs using all available data from 3.6 fim to 1250 /im 
(see Table [J). Ourl6 /Ltm flux densities are updated from 
ISevmour etall (|2007D due to newer, improved calibra- 
tions for the observing mode used. The spectra agree 
well with the photometric data at the same wavelength. 

Both 4C-H23.56 and 6C J1908-H7220 show strong con- 
tinuum luminosities consistent with the very high, 
monochromatic rest-frame 5 /im luminosities, is/xm ~ 
7 X lO^^Ltn 5 - 10 x 10"*^ ergs"!), presented by 
ISeymour et al.l ()2007t ). However, these two HzRGs dif- 
fer in the strength of their 6.2 — 7.7 /im PAH com- 
plex — which is clearly seen in 6C J1908-I-7720, de- 
spite the presence of a known powerful AGN, but is 
undetected in 4C -1-23.56. PAH emission in a power- 
ful radio galaxy is unusual; lower redshift radio galaxies 
and othe r AGN generally have much weaker PAH emis- 
sion ( e.g. iHaas et al.ll2005tlOgle et~al|[200^ [Clearv et all 



120071: IShi et al.ll2007h . In fact, the integrated PAH lu- 
minosity of 6C J1908-I-7720 is the highest of any source 
we are aware of, though the strong mid-infrared contin- 
uum conspires to make the observed 6.2 /im PAH feature 
equivalent width rather modest. 

The mid-infrared spectrum of 4C -1-23.56 is quite dif- 
ferent. Rather than showing PAH emission features, the 
observations of 4C -1-23.56 extend longward enough to 
show a flattening at rest-frame wavelengths greater than 

8.3 /im which we interpret as the onset of the 9.7 /im 
silicate absorption feature. The steep slope of the con- 
tinuum from 24 to 70 /xm. Si, ^ v^^, for both sources is 
confirmed by the high measured flux densities at 70 /im. 
The full IR SED fitting, described in the following sec- 
tion, finds a silicate absorption of T9.7^ni=0.3 ± 0.05 
for 4C -1-23.56. This value is high compared to most 
lower redshift radio galaxies of com parable luminosity 
(|Ogle et al.ll2006t IClearv et aLl|2007D. b ut low compared 
to the sources in ISaiina et al. ( 2007bD which have far 
lower radio luminosities, but reside at similar redshifts. 
The redshift of 6C J1908-f 7220 makes it impossible to 
derive a measure of the silicate absorption feature from 
Spitzer/IRS spectroscopy. 

3. FULL INFRARED SED FITTING 

Here we supplement our Spitzer observations which 
cover the 3.6 — 160 /im range with observations from 
SCUBA, CSO and IRAM which provide detections and 
upper limits for the 350 — 1250 /im range. These data 
are presented in Table [T] and the references are given 
therein. W e use the same mod el to fit these data as 
that used bv lSajina et all HoOli) and lLacv erall (|2007l ) 
which includes five components: an old (5 Gyr) stellar 
component, three dust components of different temper- 
atures (hot, cold and warm) and a PAH model. The 
age of the stellar population is unimportant as the wave- 

^ These codes are publicly available from the SSC website : 
|http : //ssc . spitzer . caltech. edu/postbcd/ir s_reductlon.html [ 



length range used in this fitting (the near-IR) varies very 
little for populations older than 1 Gyr. The far-infrared 
(far-IR) emission was fit as a modified cold blackbody 
with temperature 50 K. A warm (small grain) component 
was included with a power-law index of 2 and cutoffs 
at high and low frequency. The far-IR starburst lumi- 
nosity, LpiR, was calculated simply from the cold dust 
component. The hot dust component, associated with 
Very Small Grains, was based on a power-law with vari- 
able spectral index and was exponentially cut off at high 
frequencies with a variable cutoff frequency. The AGN 
luminosity, Lagn , was calculated from summing the hot 
and warm components. The hot dust component is then 
reddened by a dust scree n assumed to hav e the G alactic 
center extinction law of IChiar fc TielensI ()2006[ ) , likely 
to be a good approxi mation to the extrem e density en- 
vironments of AGNs (jSaiina et al.ll20d7al ). Geometries 
other than a simple screen may be less efficient at ob- 
scuring, hence this Lfir estimate is a likely lower limit. 
We then use a PAH model template based on observa- 
tions of NGC 7714 which we find to h ave no significant 
difference to M82 (jSaiina et al.ll2007aD . During the fit- 
ting several components were fixed: the amplitude of the 
stellar component was taken from Seymour et al. (2007), 
and the short-wavelength cut-off of the hot component 
and the slope of the warm component were also fixed. 
For 4C +23.56 the cold dust component was fixed to the 
maximum allowed value and for 6C J1908-I-7220 the red- 
dening was fixed to zero as we could not constrain it via 
the silicate absorption feature. 

The principle difference in the observed SEDs, pre- 
sented in Figure H is that 6C J1908-)-7220, with it's 
strong PAH emission features, is detected at both 
350 /im from CSO/SHARCII and at 850 /im from 
JCMT/SCUBA in contrast to 4C -1-23.56 (which has not 
been observed at 350 /im or detected in SCUBA obser- 
vations). While the non-detection of 4C -1-23.45 with 
SCUBA may be due to its lower redshift (and the neg- 
ative k-correction at sub- mm wavelengths), the total 
far-IR luminosity is constrained to be significantly less 
than that of 6C J1908+7220. The results of our fitting 
for both HzRGs arc presented in Tabic [2l including to- 
tal far-IR luminosities, inferred AGN luminosities, PAH 
strengths, dust mass, optical depth, and inferred SFR. 
For 4C -1-23.56 many of these quantities are upper limits. 

In addition to the detections at 24 and 850 /im, the 
upper limits at 450 and 1000 /im provide strong con- 
straints on the amount and temperature of cold dust in 
6C J1908-h7220. We obtain a starburst far-IR luminosity 
of 2.3 ± 0.4 X 10" Lq (~ 4000 ± 7OOM0yr-i, Kennicutt 
1998) where the modified black-body cold dust parame- 
ters are constrained to be Tg=50 ifc 5K and /3=1. 5 ± 0.2 
following the approach of iMcMahon et al.l (|1999f ). This 
result implies ~ 8 ± 2 x 10^ Mq of cold dust. We find an 
approximate upper-limit to the cold d ust in 4C -1-23.56 o f 
< 1.6 X 10^ Mq for Tc=50 K. Using the lPope et all (pOOTl ) 
relations, the SFR in 6C J1908+7220 inferred from our 
6.2 /im PAH luminosity is ~ 8000 M0yr^^, around a fac- 
tor of two higher than that inferred from the far-IR lu- 
minosity, albeit with considerable uncertainty due to the 
scatter in this relation. We are able to provide upper 
limits to the 6.2 /im PAH and the far-IR (starburst) lu- 
minosity for 4C-t-23.56. These limits are relatively high, 
hence the SFR in 4C +23.56 could still be at a level con- 
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sidered extrerae in the local Universe, ~ 1000 Moyr . 

4. INTERPRETATION AND CONCLUSIONS 

The clear presence of PAHs in 6C J1908+7220, a 
known powerful AGN, is consistent with results at lower 
redshift wher e PAH ernission is seen in AGN and quasar 
host spectra (|Shi et al.ll2007f) . PAHs can survive the high 
radiation environment of the AGN if shielded by an ob- 
scuring torus or located far enough from the AGN that 
they are simply stimulated and not destroyed. An HST 
WFPC2 image of 6C J1908-h7220 taken with the F702W 
filter (PI G. Miley) indicates extended, irregular star for- 
mation in the rest-frame far UV supporting the view of 
star formation occurring away from the nucleus in an ir- 
regular morphology (see Figure [3]) . The absence of PAH 
emission from 4C +23.56 is potentially due to destruction 
by the AGN, but could also simply be due to the rela- 
tively lower intrinsic SFR. No HST imaging of sufficient 
quality is available for thi s source, but it does poss ess a 
relatively large Lya halo ("Villar -Martm et al.ll2007^ im- 
plying at least a reasonably large reservoir of material 
for star formation. 

Bo t h HzRGs have deep UV spectropolarimetry (jPevI 
[19991 : iVernet et al.l[20(ll . 4C +23.56 is highly polar- 
ized (P=15.3 ± 2%), while 6C J19G8+7220 is unpolar- 
ized (P<5%) and shows broad absorption lines. The 
contribution of a scattered AGN is therefore expected to 
be much higher in 4C +23.56 than in 6C J1908+7220. 
As 4C +23.56 is almost two magnitudes fainter than 
6CJ1908+7220 at optical wavelengths (i?=23.2 for 
4C +23.56 compared to R=21A for 6C J1908 +7220 - 
[Chambers et all 119961: iPentericci eFan |2000[) . the re- 
maining rest-frame continuum emission (i.e. UV radi- 
ation from hot stars), after correction for the scattered 
AGN emission, is at least an order of magnitude fainter in 
4C +23.56 than in 6C J1908+7220. While the determina- 
tion of SFRs from the rest-frame UV continuum remains 
very uncertain in these objects, this large difference does 
suggest a much higher SFR in 6C J1908+7220, consistent 
with our IRS results. 
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The mid-IR luminosities and 70/24 /im flux density ra- 
tios (roughly equivalent to 15/7 /im rest-frame) of these 
two HzRGs are at the high e nd of the d istrib ution 
of the radio ga l axies studied by iQgle et all (|200^ and 
ISevmour et al.l (|2007l ). However, there is a marked dif- 
ference in their 8 to 24 /zm slope, beyond the fact they 
sample slightly different rest frame wavelengths. The 
steeper 8 to 24 /im slope of 4C +23.56 indicates a much 
more powerful AGN relative to the host galaxy rest- 
frame near-IR luminosity (and hence stellar mas) com- 
pared to 6CJ1908+7220. 

The mid-IR luminosities and slopes of these sources are 
also similar to high redshift 2 4 /xm selected ULIRGs stud- 
ied b y several groups {e.g., lYan et ahl 120071 : iDev et all 
l2008f ). This similarity is not surprising given the selec- 
tion of those samples (very red from i?— band to 24 ^m 
and S'24um ^ 1 mJy) , but does support our earlier com- 
ment that HzRG are representative of a larger popu- 
lation which experience brief periods of extreme AGN 
activity and obscured star formation. In fact, some of 
these 24 /im-selected ULIRGs do host radio-loud AGN 
(jSaiina et al.l l2007b[ ). but have radio luminosities sev- 
eral orders of magnitudes less than the sources presented 
here. Additionally, the correlation between the presence 
of PAH and strong emission from cold dust has been seen 
in Spitzer selected type 2 AGN (Lacy et al. 2007; Sajina 
et al. 2008). Further IRS and far-IR observations would 
provide us with an excellent opportunity to examine how 
star formation is related to the obscuration of powerful 
AGN and their radio-loud phase. 
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TABLE 1 

Photometry used in the SED fitting. 



HzRG 


wavelength 


flux density 


reference 




/^m 






J 1908+7220 


3.6 


2UU ± 20 


1 




4.5 


22y ± 23 


1 




5.8 


241 i 25 


1 




8.0 


/ion ^ A Q 
48U ± 48 


1 




12 


Q .1 A _l_ inn 
84U ± lUU 


2 




lo 


lo2U ± /U 


Q 
O 




24 


1910 + 100 


1 




70 


16200 + 1900 


1 




160 


< 63300 


1 




350 


90000 ± 15000 


4 




450 


51000 


5 




850 


13500 


6 




1250 


3000 


6 


4C +23.56 


3.6 


61 + 6 


1 




4.5 


86 + 9 


1 




5.8 


127 + 13 


1 




8.0 


424 + 043 


1 




16 


1610 + 130 


2 




24 


4390 + 200 


1 




70 


30300 + 3000 


1 




160 


< 70500 


1 




450 


< 51000 


7 




850 


< 2940 


7 
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TABLE 2 

Measured Global Properties of the HzRGs From the IR SED Fits 



HzRG 


z 


isOOMHz 


Lfir 


^AGN 


i6.2/jm 


EW6.2/jm 


T9.7/im 




SFRJ./H 






ergs~^Hz~^ 


Lq 


Lq 


Lq 


/.tm 




Mq 




6C J1908+7220 
40+23.56 


3.53 
2.48 


7.9 X 10^'' 
1.3 X 10^9 


2.3 ±0.4 X 10" 
< 6.2 X IQis 


6, 3 ±1.9 X 10" 
< 3.8 X 10^3 


8.3 ±2.0 X 10" 
< 1.1 X IQio 


0.04 ±0.01 


0.3 ± 0.05 


8 ± 2 X 10* 
< 1.6 X 10* 


4000 ± 700 
< 100 



^ The 9.7 /.im silicate feature is not covered by current observations. 
The 6.2 /zm PAH feature is not detected in this source. 
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Fig. 1. — IRS spectra of 6C J1908+7220 (left) and 4C +23.56 (right) marked by thick solid line and 1 a uncertainties indicated by dashed 
red line. Spitzer photometry is also plotted along with the model (thin line) fit to the ~ 1 — 1000 /im SED shown in Figure [2] The 
vertical arrows indicate the expected location of the strongest PAH features in the wavelength range covered by the IRS spectra (3.3, 6.2 
and 7.7^m). Note that 60 J1908+7220 exhibits PAH emission, while 40+23.56 appears to show some 9.7/im absorption, but exhibits no 
detectable PAH emission. 
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Fig. 2.— Total (1-1000 Atm) IR SEDs of 4C +23.56 and 6C J1908+220 fitted to available data. The fitted data (filled solid points and 
thick lines) are described in the text. The open points over the same redshift range as the spectra represent data not used in the fitting. 
The model (thin solid line) is comprised of an old stellar component (dashed line) and three dust components (dotted lines) of varying 
temperatures. In addition, a PAH template is included and the Galactic center extinction curve has been applied. 
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Fig. 3.— Image of 60 J1908+7220 obtained with HST/WFPC2 in the F702W band. Note the extended, diffuse emission to the NE 
at these rest-frame far-UV wavelengths, likely associated with extended star formation. The IRS observations with a wide slit completely 
cover the extended emission. 



